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Abstract 

Even if non-adiabatic perturbations are generated in multi-field inflation, 
the perturbations will become adiabatic if the universe after inflation enters 
an era of local thermal equilibrium, with no non-zero conserved quantities, 
and will remain adiabatic as long as the wavelength is outside the hori- 
zon, even when local thermal equilibrium no longer applies. Small initial 
non-adiabatic perturbations associated with imperfect local thermal equilib- 
rium remain small when baryons are created from out-of-equilibrium decay 
of massive particles, or when dark matter particles go out of local thermal 
equilibrium. 
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I. INTRODUCTION 



A recent paper [1] has shown, contrary to what is sometimes thought [2], 
that if the cosmological perturbations produced during inflation are in the 
adiabatic mode after they leave the horizon, then they stay in that mode 
through reheating and until horizon re-entry. This paper considers a comple- 
mentary question: if the cosmological perturbations produced during multi- 
field inflation are not in the adiabatic mode after they leave the horizon, then 
must they remain non-adiabatic thereafter? 

As in [1], our approach to this problem is based on the theorem that 
in all cosmological models, among the various solutions of the differential 
equations that govern perturbations in Newtonian gauge, there are always 
two independent solutions that are adiabatic throughout the whole period 
when the perturbations are outside the horizon. [3] By "adiabatic" is usually 
meant that the perturbation SS in any four-scalar S (including temperature, 

proper number densities, scalar fields, etc.) is proportional to S, where bars 
denote unperturbed quantities, and dots denote ordinary time derivatives. 
In particular, for the energy densities and pressures of all the individual 
constituents (labeled a) of the universe, we have[4] 

-^r~ = -4— is independent of a , (1) 

Pa Pa 

from which comes the term "adiabatic." The theorem of ref. [3] provides the 
further information that for all four-scalars S, the adiabatic solution has 

7-^ (2) 



and in particular 



where 



^ = ^ = -J, (3) 

Pa Pa 

lit) = f a(t)dt . (4) 



a(t) Jt 

Here a(t) is the Robertson- Walker scale factor, and ( and T are constants, 
the same for all four-scalars. (The freedom to choose T as well as ( means 
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that these are two independent solutions.) Also, according to this theorem, 
the three-scalar gravitational potentials (defined by the Newtonian gauge line 
element ds 2 = -(1 + 2$)dt 2 + a 2 (l - 2^)rfx 2 ) are [5] 

$ = = -( + HI , (5) 

(where H(t) = a(t)/a(t)), and in this mode there is no anisotropic inertia. 
(Eq. (5) can be written as the formula by which ( is usually defined[6], 
( = — ^ + 5p/3(p + p).) Observation of the cosmic microwave background 
can tell us whether the cosmological perturbations are actually in this mode 
just before they re-enter the horizon, so it is of great interest to ask, in 
what cosmological models we would expect the perturbations to be purely 
adiabatic outside the horizon. 

Before addressing this question, let's recall how current observations are 
used to check the adiabatic nature of the perturbations. It appears that dur- 
ing the whole period when the cosmic temperature is in the range 10 9 °K 
T 3> 3000 °K, the constituents of the universe consist of a non-relativistic 
nucleon-electron plasma, plus photons, plus highly relativistic free neutri- 
nos, plus non-relativistic cold dark matter. To a good approximation, energy 
is separately conserved in each of these fluids, so each satisfies the energy 
conservation conditions 

p- a + 3H(p a +p a ) = (6) 

and, outside the horizon, 

5p a + 3H{5p a + 5 Pa ) = 3{p a + p a )V . (7) 

where a runs over the four fluids — cold dark matter, neutrinos, photons, 
and the electron-nucleon plasma. For each of these fluids the pressure p a 
is either negligible or a function only of the energy density p a , so it follows 
from Eqs. (6) and (7) that [7] 

Spa * + (a , x 

Pa H 

where the ( a are various constants. (There are actually three separate con- 
stants for the three neutrino flavors.) According to Eq. (1), an adiabatic 
mode would have all ( a equal. In fact, according to Eqs. (3) and (5), they 
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would all be equal to (. The constants ( a affect the further evolution of the 
constituents of the universe after horizon re-entry and up to the time of last 
scattering, and so can be measured by observations of anisotropies in the cos- 
mic microwave background. So far, there is no evidence for unequal ( a , and 
some indications that they are equal for cold dark matter and photons. [8] 

The most usually considered case in which we would expect purely adia- 
batic perturbations is inflation with a single scalar "inflaton" field. It is well 
known that in this case the cosmic perturbations are purely adiabatic during 
inflation after horizon exit (whether or not the "slow roll" conditions are 
satisfied), and reference [1] shows that they then remain adiabatic through 
the period of reheating and as long as the perturbations remain outside the 
horizon. In this case, all ( a would be equal. With several scalar fields the 
perturbations are not in general purely adiabatic during inflation, and it 
seems to be widely expected that this would show up in unequal values for 
the constants ( a . 

But there is a plausible scenario in which non-adiabatic fluctuations pro- 
duced in multi-field inflation become adiabatic after reheating, and then re- 
main adiabatic until horizon re-entry. Suppose that reheating leaves the 
universe in a state of local thermal (including chemical) equilibrium, with no 
non-zero conserved quantities. This is what is generally assumed in conven- 
tional theories of cosmological baryogenesis. In this case the perturbations 
are described by just three degrees of freedom: the temperature fluctuation 
5T, the scalar potential $ = and the velocity potential fluctuation 8u. 
As described in [3], these variables are governed by three first-order differ- 
ential equations, so they must have three independent solutions, but there 
is also a constraint relating 5T, and Su, leaving just two independent 
solutions. Since there are always two independent adiabatic solutions, the 
solutions here are all adiabatic. (This argument would not apply if one or 
more scalars remained decoupled from the matter in thermal equilibrium [9] , 
as in curvaton models [10], theories of scalar field baryogenesis [11], and in 
some theories of cosmological moduli [12] or axions [13].) 

Later the universe develops a non-zero baryon number; neutrinos and cold 
dark matter decouple from other particles; and local thermal equilibrium is 
lost. The arguments of [3] then no longer apply, which is why the field 
equations in the era from e + e _ annihilation to last scattering have solutions 
(8) with different ( a . Nevertheless, as long as the wavelength remains outside 
the horizon the adiabatic solution (2)-(5) is always a solution, and since 
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the cosmological perturbations were described by this solution during the 
assumed early era of local thermal equilibrium, they remain purely adiabatic. 

Just as in [1], there is a weak point in this argument. Local thermal equi- 
librium is never perfect, so there are always other degrees of freedom besides 
ST, and Su, and hence other solutions of the field equations besides the 
adiabatic solutions. We shall show in Section II that any non-adiabatic con- 
tributions to cosmological fluctuations will decay as the universe approaches 
a state of local thermal equilibrium with zero chemical potentials. But even 
if these non-adiabatic solutions make only a tiny contribution to cosmological 
perturbations as this state is approached, how do we know that they do not 
grow rapidly again as local thermal equilibrium is subsequently lost, as for 
instance by the decoupling of dark matter particles, or as non-zero chemical 
potentials appear in cosmological baryosynthesis? 

It seems that this question must be addressed on a case by case basis. 
Here we will consider it in two familiar contexts: the survival of cold dark 
matter particles as they go out of local thermal equilibrium, treated in Section 
III , and the appearance of a non-zero chemical potential through the out- 
of-equilibrium decay of a massive exotic particle, considered in Section IV. 

II. APPROACH TO ADIABATICITY 

Suppose that after inflation the dominant constituent of the universe is 
a heat bath in local thermal equilibrium with no chemical potentials, but 
that, as a remnant of multi-field inflation, there is also present a species of 
particles with number density n that is not even approximately described 
by the adiabatic condition (2). The rate of change of this number density 
in co-moving inertial frames will be some function Y{n,T) of n and of the 
temperature T of the heat bath 2 

(n«% = y. (9) 

2 Here u M is the velocity four- vector, normalized so that g^u^u" = — 1, with unper- 
turbed components u° = +1 and u % = 0. Because outside the horizon spatial gradients 
can be neglected, it is unnecessary to state whether is the velocity four-vector of the 
particles with number density n, or of the heat bath energy-momentum tensor, or of the 
total energy- momentum tensor — all we need to know about 5u^ is that Su° = 5uq = — 
which follows from the normalization condition on 
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To zeroth order in perturbations this is 



n + 3Hn = Y (10) 

while to first order 

5n + 3H5n - 3n^ = 5Y + $Y (11) 

For simplicity, we will assume that there were so many particle species in 
equilibrium that n can be ignored in the evolution of \1/ and T. According to 
the argument quoted in Section I, it follows that 5T and \1/ are given by the 
adiabatic conditions: 

= -X , $ = ^ = Hl-C ■ (12) 

T 

It follows that the perturbation to the rate of change of n takes the form 

dY 

5Y + Y1= —(Sn + nX) , (13) 
on 

the terms involving dY/dT canceling. Adding the time derivative of X times 
Eq. (10) to Eq. (11) and using Eqs. (12) and (13) gives 

±(6n + Xh) = ^-3H+^j(Sn + Xh). (14) 

This must be compared with the rate of decrease of n itself, given by Eq. (10). 
From Eqs. (10) and (14) we have 



dt 



Sn + nl 



n 



dY Y 

on n 



The right-hand side may have any sign, but we know that when n is near the 
value nEQ(T) that it would have in local thermal equilibrium at a temperature 
T, the rate of change of n due to all causes apart from the expansion of the 
universe — decay, annihilation, etc. — must take the form 

Y - -A(T) (n - n EQ (T)) (16) 
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where A is some rate that must be positive in order that n in Eq. (10) should 
be attracted toward tzeq when the expansion of the universe can be neglected. 
For A(T) 3> H the expansion of the universe can be neglected, and so as the 
density n approaches the value it would have in thermal equilibrium, the 
right-hand side of Eq. (15) approaches the negative value — A(T), and the 
ratio (5n + nX) / n, which can be taken as a dimensionless measure of the 
departure from adiabaticity, decays exponentially with rate A(T). We will 
see in Sections III and IV that this ratio is the quantity that sets the scale 
for possible inequalities in the observable quantities ( a . 

We have seen that the non-adiabaticities associated with multi-field in- 
flation become small if inflation is followed by a period of local thermal equi- 
librium with no non-zero chemical potentials. Now we have to see if they can 
be revived when the universe subsequently goes out of thermal equilibrium. 

III. DARK MATTER 

Next we consider the possible non-adiabatic contributions to cosmological 
perturbations that might be caused by the departure from local thermal 
equilibrium that occurs when the annihilation rate of heavy cold dark matter 
particles can no longer keep up with the rate that these particles would 
disappear in thermal equilibrium. [14] 

We assume here that all matter and radiation except the cold dark matter 
particles are in local thermal equilibrium with each other, and that this 
thermal bath dominates the energy density of the universe. Then the density 
n of the cold-dark matter particles will again be governed by Eq. (15). As 
we have seen in Section II, if there is a sufficiently long period when the 
rate Y characterizing interactions of the cold dark matter particles with the 
thermal bath is much larger than the expansion rate, then by some time 
t\ the ratio \Sn + n1\/n will have become exponentially small. Eventually, 
as the temperature falls below the dark matter particle mass, the fractional 
rate of decrease of tzeq will become greater than Y, and n will consequently 
become greater than ueq, eventually leaving over a remnant of dark matter 
that has survived to the present. The question is whether the small non- 
adiabatic perturbation at time t\ is amplified while the dark matter goes out 
of equilibrium. 
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Even when the cold dark matter density n begins to differ appreciably 
from its equilibrium value ueq, Eq. (15) will still apply but Y will no longer 
be given by Eq. (16). If we assume that dark matter particles are annihilated 
only in binary collisions, and created only in pairs, then Y will be of the form 



Y(n,T) = -R(T)(n 2 -n 2 EQ (T)) , 



(17) 



where R is a positive rate constant, equal to the average over the dark matter 
velocity distribution of the product of the dark matter particle velocity and 
the annihilation cross section. Then Eq. (15) becomes 



dt 



Sn + nX 



n 



-R(f) 



' n 2 + n 2 



EQ 



cry 



n 



(18) 



At late times n 3> tieqCO, so the right-hand side approaches —R(f)n. The 
rate constant R(T) approaches a constant at low temperature, and n asymp- 
totically goes as a -3 , so the time-integral of the right-hand side of Eq. (18) 
converges for large t. Therefore at late times 

Sn + hX] \ f°° (n 2 +n EQ (fy 



Sn + nX 



n 



n 



cxp 



ti 



R(T) 



dt 



n 



(19) 



We see that the departure from adiabaticity at late times is less than at 
the end of the period when the cold dark matter particles were in thermal 
equilibrium with everything else, though only by a finite factor. 

At late times n — > —3Hn and HX = ( + \1>, so Eq. (19) may be written 
as in Eq. (8): 

5p D 5n f( D + ^\ 



Pd 



n 



H 



with ( D now given by 

Sn + nX 



c 



D 



n 



cxp 



J i 



ti 



R(f) 



n 2 + n 2 EQ (TY 
fi 



dt 



(20) 



(21) 



As we have seen, the quantity [(Sn + nX)/n} 1 is exponentially small if there 
is a sufficiently long period when the rate Y characterizing interactions of 
the cold dark matter particles with the thermal bath is much larger than the 
expansion rate, and Eq. (21) shows that then the departure of the observ- 
able quantity (o from the value ( that it would have for a purely adiabatic 
perturbation is even smaller. 
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IV. BARYOGENESIS 



Here we will consider the simple original model of cosmological baryoge- 
nesis[15], with the later modification [16], that the out-of equilibrium decay 
of some heavy exotic particle X produces a non-zero value of the quantity 
B — L, which is subsequently processed by B — L-conserving non-perturbative 
effects of the electroweak interactions[17] into some definite proportions of 
B and L, depending only on the numbers of generations and scalar dou- 
blets^]. It is assumed that at some time t\ after inflation a nearly perfect 
local thermal equilibrium has been reached, with no chemical potentials, but 
that subsequently the actual rate of disappearance of the X particle (and its 
antiparticle, if distinct) does not keep up with the rapid decrease that would 
occur in thermal equilibrium as the temperature falls below the X mass. In 
the original model [15] of direct baryosynthesis the X particles were assumed 
to have distinct antiparticles, but with equal number densities nx- In the 
case [16] of leptogenesis the X particles are usually assumed to be identical 
with their antiparticles, and their density will again be denoted nx- We 
will assume that there are so many particle species in thermal equilibrium 
that the density of the X particles makes a negligible contribution to the 
gravitational field and to the evolution of the temperature of the particles in 
equilibrium, so that the fluctuations in the temperature and metric are given 
by Eq. (12). Then the degree of non-adiabaticity of this number density will 
again be given by an equation like Eq. (15), in which Y is now the rate of 
change of nx in co-moving inertial frames. As before, this gives 



a subscript 1 denoting the time t±, when the decay of the X particles has not 
yet begun to produce any appreciable net density of B — L. But once the 
temperature falls below the X-particle mass and the disappearance rate of 
the X-particles falls below the expansion rate, the density nx will be much 
larger than its equilibrium value, and Eq. (16) will not apply. 

Suppose that the disappearance of a single X particle (together with its 
antiparticle, if distinct) produces on average a net value b of B — L. The rate 
of change of the density n B _ L of B — L in co-moving inertial frames is then 




(22) 




(23) 
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To zeroth order in perturbations, this gives 

h B - L + 3Hn B - L = -bY , (24) 

while to first order 

biiB—L + 3H5n B ^ L - 3n B ^ = -b{SY + $?) . (25) 
The sum of Eq. (25) and the time derivative of X times Eq. (24) is 

(cm B -L + ^b-lI) + 3H(8tib-l + ^b-lI) = -b -^—(Sn x + 1h x ) • (26) 



dt 

The solution is 

8n B -L + fiB-iZ = b 



exp 



* dY 



i dn 



dt 



x 



Sn x + n x I 



(27) 



a subscript 1 again denoting quantities evaluated at time t±. Also, the solu- 
tion of Eq. (24) is 



n B - L = ~b 



n x - 



(28) 



Thus the departure from adiabaticity in the B — L density is 



5n 



B-L 



+ lh 



B-L 



n B - L 



(Sn x + Xh x ) l 
n xl - (a/a 1 ) 3 n Xi 



cxp 



* dY 



dn 



dt 



x 



(29) 



Asymptotically only decay contributes to the disappearance rate Y, so Y — > 
—Fxnx (where T x is the decay rate of the A-particle), and therefore the 
exponential in Eq. (29) becomes negligible compared with one and (a/ai) 3 n x 
becomes negligible compared with n X i- Also, after the temperature drops 
below the electroweak scale 300 GeV, the values of the zeroth and first- 
order terms in the baryon density n B and 5n B are proportional respectively 
to n B ^L and 5n B _ B with the same constant factor, so 



Sn B + n B l 8n B _ L +lh B ^ L 



n B 



n B - L 



'5n x + Xn x 
n x 



(30) 
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Equivalently, since at late times n B oc a 3 and HT + 

5n B /^ + Cb-l\ 



n B 



H 



(31) 



where (b is the constant 



Cb = C + 



5n 



x 



n x l 



(32) 



As we saw in Section II, if there is a sufficiently long period ending at time 
t\ during which the X particles are decaying and being recreated by the 
heat bath of other particles at a rate much greater than H, then the second 
term on the right-hand side of Eq. (32) will be exponentially small. Hence, 
on these assumptions, the observed constant (b will be very close to the 
common value ( of the other ( a . 



V. CONCLUSION 



If we assume that the results of Sections III and IV are typical of the 
general case, then the existence of a post-inflation period of nearly perfect 
local thermal equilibrium would rule out any appreciable present departures 
from a purely adiabatic perturbation. Thus the observation of a departure 
from a purely adiabatic perturbation, such as a measurement of different 
values for the ( a in Eq. (8), would be evidence not only that non- adiabatic 
perturbations are generated in inflation, but also that there was no era of per- 
fect or nearly perfect local thermal equilibrium when all conserved quantum 
numbers vanished. 
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